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Abstract 

The modest addition of the dimension-5 term A {Hu-H^)'^ /M to the superpotential of the minimal 
supersymmetric standard model (MSSM) originated from physics beyond the MSSM (BMSSM) has 
a significant impact on the scenario of the Higgsino-dominated neutralino state being the lightest 
supersymmetric particle (LSP). It increases the mass difference between the LSP and the lighter 
chargino as well as that between the LSP and the second-lightest neutralino. This enhances the 
LHC discovery potential of the chargino and neutralino decays, producing more energetic charged 
leptons or pions than the decays without the BMSSM corrections. Furthermore, the coannihilation 
between the lighter chargino or second-lightest neutralino and the LSP is reduced substantially such 
that the LSP mass does not have to be very heavy. Consequently, an almost pure Higgsino LSP 
with its mass ~ 100 GeV in the BMSSM can account for all the relic density of cold dark matter 
in the Universe unless tan (3 is too large. 
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Introduction - The presence of cold dark matter (CDM) in our Universe is now well 
established by the very precise measurement of the cosmic microwave background radiation 
in the Wilkinson Microwave Anisotropy Probe (WMAP) experiment [l]. A nominal 3a 
range of the CDM relic density is 

^CDM = 0.105 toxm 5 (1) 

where /i is the Hubble constant in units of 100 km/Mpc/s. 

One of the most appealing and natural CDM particle candidates is provided by super- 
symmetric models with i?-parity conservation This -R-parity conservation ensures the 
stability of the lightest supersymmetric particle (LSP) so that the LSP can be CDM. The 
LSP is in general the lightest neutralino, a linear combination of neutral electroweak (EW) 
gauginos and Higgsinos. Since the LSP nature depends on its compositions, its detection 
can vary a lot. 

An interesting scenario is the Higgsino-like LSP, which can arise from a number of su- 
persymmetry breaking models, e.g., focus-point supersymmetry models 3| or whenever the 
yU parameter is much smaller than the Bino and Wino masses ^ . In the minimal supersym- 
metric standard model (MSSM) 5|], the Higgsino-LSP scenario implies nearly-degenerate 
Higgsino states: Coannihilation is too efficient so that the observed CDM relic density 
requires a rather heavy Higgsino state with mass around 1-1.2 TeV 61]. Moreover, the 
mass degeneracies between the LSP (xi) ^-iid the lighter chargino/second-lightest neutralino 
{X1/X2) generate too soft decay products for detecting the states xf X2 ^^e LHC. 
Thus, the Higgsino-LSP in the MSSM posts a difficult scenario at the LHC. 

In this Letter we show that the modest addition of a dimension-5 term A ■ HctY /M to 
the MSSM superpotential, a scenario beyond MSSM (BMSSM) alleviates the difficulties 
of the Higgsino-LSP scenario. It has been discussed that non-renormalizable and high- 
dimensional operators in new physics models can yield important consequences in low energy 
phenomenology As shall be demonstrated in the following, this dimension-5 term. 



A {Hu-Hdf/M, lifts up the de generacy between the states x? and Xi /X2 thus enhancing 
the discovery potential of the chargino and neutralino decays with more energetic charged 
leptons or pions than the decays in the MSSM. In addition, the coannihilation of the LSP 
with xt oi' X2 is reduced substantially such that a light Higgsino-LSP with a mass around 
100 GeV can accommodate the WMAP data on the CDM relic density. 
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BMSSM -Albeit many virtues of low energy supersymmetry (SUSY), fine-tuning in 
the lightest Higgs boson mass my, motivates additional degrees of freedom to the MSSM 
[lo| . New interactions beyond the MSSM at the TeV scale M may be encoded in higher- 
dimensional operators. Recently, Dine et al. 7|] have shown that the most general dimension- 
5 superpotential term for the MSSM Higgs sector is 

= 4 {Hu-hX , (2) 



M 

with the SU(2) contraction ■ = H^H^ — H^H^ for the up-type and down-type Higgs 
doublet superfields, Hu and Hd, respectively. This dimension-5 operator has been shown to 
raise easily the lightest Higgs boson mass above the LEP bound without loss of naturalness 



[III- 



Another dimension-5 operator, which breaks SUSY and affects the Higgs spectrum, is 

SdZ^{H^.H,)\ (3) 

where Z = 6'^msusY is the spurion field with the SUSY breaking scale msusY 0] • If "^susy — 
the correction to mh comes dominantly from the supersymmetric operator in Eq. ([2]) 
rather than that in Eq. ([3]). The correction is given to leading order in the dimensionless 
parameter e = Xfi/M hj 



5m\ = ev^ 



1 + ISifi + ^ ^ 



(m\ — m|)2 -I- ^m\m\s\p 



(4) 



m\ — m\ tan/3 ytan^/3 

where v 246 GeV is the Higgs vacuum expectation value and the second expression holds 
for large tan/3. For simplicity we take the CP-conserving framework and decoupling limit 
of mj^ ^ mz in the following. Note that the BMSSM correction is inversely proportional to 
tan/3 in contrast to the conventional radiative corrections. The correction bnih normalized 
by vtih around the LEP bound is roughly bmhjvtih ~ l^ej tan/3 for rtiA ^ vtiz- For e = 0.05 
the correction can be as large as 50% (10%) for tan/3 = 2 (10). Therefore, the LEP bound 
is easily satisfied by a positive e ~ 0.1. 

The interaction terms that involve only the Higgsino fields [Hu^d) and Higgs fields {Hu^d) 
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are given by 



-2 ( Hu ■ Hd) ( Hu ■ Hd) + ( Hu ■ Hd) + (Hu ■ Hd 



E.G.. 



(5) 



After EW symmetry and SUSY breaking, the modified neutralino mass matrix A^at in the 
{B, W\ H% Hi) basis reads: 



M 



N 



/ Ml 
M2 

-mzswcp mzcwcp 



-mzSwCfs 
■mzcwcp 



mzSwS/3 
-nizcwsp 



\ 



\ mzswsp -mzCwSfs -ji + jjv^cpsp c^p 



(6) 



M 



c 



(7) 



and the modified chargino mass matrix Aic iii the {W , H } basis reads: 

where = sin /3, sw = sin^^yj etc. To leading order in e, the BMSSM effects on the masses 
of x?, ^iid Xi are, in the hght Higgsino case (Mi^2 ^ 12\ 



Xl,2 



m~± ~ 

Xi 



m 



2/3 



2^' 



2^2/3^ 



- sign(/i) 



sign(/i) 



2M1M2 



where in the first equation the upper (lower) sign is for x\ (X2) mass. With increasing e both 
the Xi aiid Xi masses decrease; the X2 mass increases for large tan (3 > A but decreases for 
small tan/5 < 4, if the mv//Mi^2 correction is ignored. However, since the LSP mass drops 
faster than the Xi mass with increasing e, sizable mass differences between Xi a^id x\ as 
well as between X2 and x^ are developed. Due to the decreasing mass of xt with e, the lower 
mass bound of m~± > 94 GeV [3] can constrain the BMSSM light Higgsino-LSP scenario. 
If yU is negative, the second term in the expressions of m~± in Eq. ([8]) slows down the lighter 
chargino mass, which leads to larger mass difference between Xa and xf . Therefore, the 
negative yU. case accommodates larger parameter space to explain all the WMAP data by 
the Higgsino-LSP. However we note that the combined analysis for the anomalous magnetic 
moment of the muon as well as 6 — 57 prefers a positive /i [l^ . In what follows, therefore, 
we take the case of positive /i. 



Collider detection - It is well known 12|] that in the Higgsino-LSP scenario the mass 
degeneracy among x5, ^ind xf renders their detection at colliders extremely difficult, even 
though radiative corrections can increase the mass difference by a few GeV |15|^. The decay 
products of Xi' or X2 ^'^'^ soft for their detection. In the BMSSM, the new contribution 
from the dimension-5 term in Eq. ([2]) can alleviate the Higgsino-LSP detection problem 
significantly by inducing sizable mass splittings between Xi and xf /X2- ^^e other hand, 
the Higgsino fraction of the LSP remains high enough to call the LSP a pure Higgsino as long 
as the gaugino masses are large. The effect of e on the neutralino mixing matrix A^, which 
diagonalizes the neutralino mass matrix as N*Ai]yN'^ = diag(m^o, . . . , m^o), corresponds 
to the rotation of and components. Therefore the LSP Higgsino fraction, Pg^ = 
l^iaP + |^i4p; remains intact by the change of the e parameter. In the limit of large 



gaugino masses, we have almost pure Higgsino LSP since Pg 
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Figure 1: The mass difference Amo = m~o — m-o (solid) between the second lightest neutralino 
and the LSP and the mass difference Am± = rn~± — m~o (dashed) between the lighter chargino 
and the LSP as a function of the BMSSM correction parameter e. 



^ The radiative corrections to the mass sphttings can be as large as 10 GeV if the stop mixing angle and 
mass splitting are large. See Ref. 16 1. 
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In Fig.d], we show the mass splittings, Amo = 1^X2 ^"^x? Am± = m~± — m^o, due to 
the BMSSM corrections as a function of e for a specific choice of parameters: M2 = 2Mi = 1 
TeV, /i = 120 GeV and tan/3 = 3. In this parameter set, the LSP mass is 70 — 110 GeV. 
For e = 0.05-0.1, the mass difference between xf cind Xi is about 15 — 30 GeV and the 
mass difference between X2 Xi is about 15-40 GeV. These mass splittings are much 
larger than those due to radiative corrections, which are typically a few GeV. Such sizable 
mass differences can help us detect the lighter chargino or second-lightest neutralino by 
tagging more energetic charged leptons or jets in the decay products. Consequently, the 
phenomenological impact of the BMSSM corrections on the SUSY search at the LHC is 
expected significant. Nevertheless, we do not perform any full-fledged analysis in the present 
work, expecting that such a comprehensive analysis will lead to almost the same physical 
conclusions as those described above. 

Dark Matter - In the MSSM Higgsino-LSP scenario, the strong coannihilation due to 
mass degeneracy pushes the Higgsino mass rather high, about 1-1.2 TeV, to account for 
the CDM relic density in Eq. ([1]). We note, in passing, that the MSSM radiative corrections 
affect the relic density rather mildly 17|]. In the BMSSM, however, the large mass differences 
of Amo and Am-t can suppress the coannihilation effectively; a much lighter Higgsino-LSP 
can account for the CDM relic density. 

Another compelling feature in the BMSSM arises from the modified Xi'Xr^ coupling 
which is proportional to (jA'^isp — |A'^i4p). The e dependence can be easily seen from the 
first row of the matrix to leading order in the BMSSM corrections, as 

Nu -(0,0, 1±|^, , (9) 



V2 ' V2 

where Eh = ew^C2/3/(4/i^), and we ignore small terms of 0{mw/Mi) as well as an overall phase 
fisl . In the MSSM {e = 0), the light Higgsino-LSP scenario implies an almost vanishing 
XrXr^ vertex. In the BMSSM, the modified neutralino mixing matrix in Eq. ([9]) leads 
to a sizable XrXrZ vertex, linearly proportional to e. Therefore, the annihilation process 
XiXi ^ Z ^ f f can be enhanced by the BMSSM corrections, giving a profound effect on 
the relic density of the LSP. 

For a simple quantitative estimate we use a useful formula for the CDM relic density 

n,h^ ^ . (10) 
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We include all the 2^2 self-annihilation and coannihilation processes in calculating the 
effective annihilation cross section acs- Since the mass difference ArriQ is substantially larger 
than Am± for e ~ 0.1 (see Fig.[T]), we ignore the X1X2 coannihilation in estimating {ccsv) and 
use the following formula for a crude estimate of the thermally-averaged effective annihilation 
cross section, taking into account the coannihilation from xf- 



We take the freeze-out temperature Tf = m^o/25 and the relative velocity Vij = 0.3 during 



XiXi ^ hh, Z^h, Z^Z^, W^W~ , ff. For the coannihilation cross section (cr-o^-) we include 
the processes XiXi ~^ W~h^, W~Z^, W~'-f, ff' [l^. In the following numerical analysis we 
fix ttLq^ = rriqj^ = mj = —A = rriA = 1 TeV and include all the radiative corrections to the 
masses of neutralinos, charginos, and the Higgs bosons. 

In Fig. [21 we show the tan/3- and Mi-dependence of the relic density Vl^h"^ versus e 
with the fixed /i = 120 GeV. In Fig. [2](a), we fix tan/9 = 3 and take three typical values 
of Ml = 300, 500, 700 GeV. As expected, the relic density increases with e due to the 
suppressed coannihilation. In addition, we observe that decreasing Mi reduces the LSP relic 
density. This is because the effective annihilation cross section o"cfT, dominated by the process 
XiXi ~^ ff for sizable e, is enhanced by the stronger Bino-Higgsino mixing for smaller Mi. 
Nevertheless the Mi-dependence of the relic density is rather mild. 

On the other hand, the tan /3-dependence of fi^/i^ is strong, as can be seen from Fig. [2](b). 
We take tan/3 = 2, 3, 4, 5, 10 for the fixed /i = 120 GeV and Mi = 500 GeV. The tan/3 = 2 
case has the curve terminated at a large e as m~± gets below the experimental bound, 
m~± > 100 GeV. Up to e ~ 0.08, the relic density is increasing with e, because of more 
suppressed coannihilation. And is larger for smaller tan/3, which can be attributed 
to the tan /3-dependence of the masses in Eq. ([8]): Small tan/3 ^ 1 maximally reduces the 
LSP mass, and thus enhances the mass difference Am±, suppressing the coannihilation. 
The slope of increasing fi^ around e ~ 0.08 is very steep, due to the kinematic closure of 
XiXi W~^W~ mode {m^o is decreasing). Another interesting feature is that after e ~ 0.08 
the relic density in the small tan/3 case turns its direction and decreases. This is because. 




3/2 



(11) 



the freeze-out. For the self-annihilation cross section (a^o^o) we consider the processes 
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Figure 2: The Higgsino-LSP relic density Q.^h? as a function of e: (a) for the fixed tan (3 = 3 with 
Ml = 300,550,700 GeV, and (h) for the fixed Mi = M2/2 = 300 GeV with tan/3 = 2,3,4,5,10. 
We set n = 120 GeV and other SUSY mass parameters to 1 TeV. 

as e increases, the Xi'Xr^ vertex becomes stronger and thus the self-annihilation process 
XiXi ^ Z ^ ff is enhanced. Finally we observed that the WMAP data of VL^h^ ~ 0.1 
with positive /i strongly prefer small tan/3. If tan/3 > 4, the WMAP data can be only 
partially explained. Of course, if is negative, a much larger parameter space is allowed for 
the WMAP data. 

In order to understand the behavior of the relic density against e (especially for the 
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Figure 3: The effective self and coannihilation cross sections. TheW~^W mode is from the XiXi 
self-annihilation while the W^Z/W^^ modes from the Xixf coannihilation. 

bumpy shape), we calculate the contribution of each channel for small tan/3 = 2 case. 
Figure [3] shows the effective annihilation cross section a^s as a function of e. Here we 
set /i = 120 GeV, and Mi = 500 GeV. We take into account the thermal suppression 
factor due to the mass difference. The effective Xixf coannihilation cross section is defined 



decreases. For e > 0.072 the W~^W~ mode is kinematically closed. The self-annihilation 
process XiXi ~^ ff is enhanced due to the XrXi-Z vertex being stronger with increasing e. 
On the other hand, the coannihilation channels, dominant for small e, become suppressed 
with larger e due to larger mass splittings. As a whole, we have the bump-shaped distribution 
of the relic density as a function of e in Fig|2]^b). One technical issue can arise when we 
calculate the thermally averaged annihilation cross section near W~^W~ thresholds. Since 
we fix the relative velocity, there is some discrepancy in the cross section which is sensitive 
to the actual relative velocities near the threshold. A more detailed analysis based on the 
exact Boltzman equation will be required for more accurate estimates of the cross sections, 
which is, however, beyond the scope of the present short report. 

Another important experimental test for the Higgsino-LSP scenario is the spin- 
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Figure 4: As functions ofe, we present the relic density Q^h? in units of 0.1, the lightest CP-even 
neutral Higgs boson mass rufi, the LSP mass, and the lighter chargino mass in units of 100 GeV, 
and the Higgsino fraction of the LSP. We fix fj, = 120 GeV, tan /3 = 3, and Mi = 500 GeV. 
Radiative corrections to the neutralino and chargino masses as well as to nih are also included. 

independent scattering cross section cr^p of the LSP with nucleons. We have calculated 
the scattering cross section based on the input parameters given in Ref. 20|. As shown 
in Fig.m if e is larger than 0.02, the elastic scattering cross section is reduced far be- 
low the current limits, mainly because the lightest Higgs boson mass rrih increases with 
the BMSSM corrections. Here we briefly comment on the spin-dependent scattering cross 
sections. The spin-dependent scattering process, which is experimentally more difficult to 
extract because of the lack of coherent enhancement unlike the spin-independent process, 
is mainly mediated by the Z boson. It is enhanced by the stronger XrXi-Z coupling with 
increasing e. 

Figure m summarizes all of our findings of the light Higgsino-LSP scenario in the BMSSM. 
For the MSSM parameters, we have chosen moderate values which can explain the WMAP 
data: fi = 120 GeV, tan (3 = 3, and Mi = 500 GeV. We show, as functions of the parameter 
e, the relic density Qg in units of 0.1, the lightest Higgs boson mass m/^, the LSP mass 
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m^o, the xf mass (all masses in units of 100 GeV), the Higgsino fraction of the LSP, 
and the spin-independent scattering cross section 0"^^ in units of 10~'^^cm^. As e increases, 
the Higgs mass increases but both m^o and m~± decrease. However, the mass difference 
Am± increases while the Higgsino fraction of the LSP stays high ~ 98%. We see a 
dramatic enhancement of the Higgsino-LSP relic density fi^ h"^ near e = 0.08 or larger, 
where XiXi ~^ W~^W~ channel is kinematically closed. For our parameter choice, the light 
Higgsino-LSP of its mass about 82 GeV, when e ~ 0.085, can explain all the observed CDM 
relic density in the Universe. For other choices of parameters we can still accommodate 
a light Higgsino-LSP as a primary candidate for the CDM relic density observed by the 
WMAP unless tan (3 is large. 

Conclusions - The fine-tuning of the lightest Higgs boson mass ruh in the MSSM motivates 
additional interactions beyond the MSSM around the TeV scale. We have checked that, in 
the effective Lagrangian approach, the least-suppressed dimension-5 operator A {Hu-HdY /M 
added to the MSSM superpotential can usually increase sufficiently for a moderate 
tan/3 and it can significantly affect the light Higgsino-LSP scenario. It lifts up the mass 
degeneracy between the LSP and the lighter chargino and that between the LSP and the 
second-lightest neutralino as much as a few tens of GeV. As a result, it is expected to improve 
significantly the chance of detecting the decay products of charginos and neutralinos at the 
LHC. Another important impact is on the Higgsino-LSP dark matter. Since the mass 
splittings suppress the coannihilation processes, the WMAP narrow band for the CDM relic 
density can accommodate a light Higgsino-LSP particle of its mass around 100 GeV for 
rather small tan (3. 
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